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Abstract 


The  Tise  of  a  blade  vortex  interaction  noise  prediction 
scheme,  based  on  CAMRAD/JA,  FPR  and  RAPP, 
quantifies  the  effects  of  errors  and  assumptions  in  the 
modeling  of  the  helicopter’s  shed  vortex  on  the  acous¬ 
tic  predictions.  CAMRAD/JA  computes  the  wake  ge¬ 
ometry  and  inflow  angles  that  are  used  in  FPR  to  solve 
for  the  aerodynamic  surface  pressures.  RAPP  uses 
these  surface  pressures  to  predict  the  acoustic  pres¬ 
sure.  Both  CAMRAD/JA  and  FPR  utilize  the  Biot- 
Savart  Law  to  determine  the  influence  of  the  vortical 
velocities  on  the  blade  loading  and  both  codes  use 
an  algebraic  vortex  model  for  the  solid  body  rotation 
of  the  vortex  core.  Large  changes  in  the  specifica¬ 
tion  of  the  vortex  core  size  do  not  change  the  inplane 
wake  geometry  calculated  by  CAMRAD/JA  and  only 
slighty  affect  the  out-of-plane  wake  geometry.  How¬ 
ever,  the  aerodynamic  surface  pressure  calculated  by 
FPR  changes  in  both  magnitude  and  character  with 
small  changes  to  the  core  size  used  by  the  FPR  calcu¬ 
lations.  This  in  turn  affects  the  acoustic  predictions. 
Shifting  the  CAMRAD/JA  wake  geometry  away  from 
the  rotor  plane  by  1/4  chord  produces  drastic  changes 
in  the  acoustic  predictions  indicating  that  the  predic¬ 
tion  of  acoustic  pressure  is  extremely  sensitive  to  the 
miss  distance  between  the  vortex  and  the  blade  and 
that  this  distance  must  be  calulated  as  accurately  as 
possible  for  acceptable  noise  predictions.  The  inclu¬ 
sion  or  exclusion  of  a  vortex  in  the  FPR-RAPP  cal¬ 
culation  allows  for  the  determination  of  the  relative 
importance  of  that  vortex  as  a  BVI  noise  source.  Be¬ 
ing  able  to  identify  which  vortex  is  responsible  for  the 
largest  noise  sources  will  lead  to  improved  BVI  noise 
reduction  methods. 


This  paper  is  declared  a  work  of  the  U.S.  Government  and 
is  not  subject  to  copyright  protection  in  the  United  States. 
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Introduction 

Under  certain  flight  conditions,  the  helicopter  main 
rotor  blade  passes  through  or  very  close  to  the  trailing 
tip  vortices  previously  generated  by  the  rotor  blades. 
This  phenomenon  is  referred  to  as  blade- vortex  inter¬ 
action  (BVI)  and  occurs  mainly  for  powered  descent, 
approach  and  nap-of-the-earth  maneuvers.  The  blade- 
vortex  interactions  produce  impulsive  changes  in  blade 
surface  pressure  distributions  which  lead  to  vibration 
and  noise  problems.  The  BVI  noise  produced  is  highly 
annoying  and  limits  the  public  acceptance  of  many 
helicopter  applications.  The  reduction  of  BVI  noise 
is  therefore  a  major  concern  to  both  the  rotorcraft 
industry  and  the  prospective  buyer. 

The  first  step  in  reducing  BVI  noise  is  to  under¬ 
stand  the  noise  generating  mechanisms.  One  of  the 
first  issues  in  the  modeling  of  helicopter  acoustics  was 
the  scalability  of  the  BVI  noise  sources.  It  was  de¬ 
termined  that  small-scale  wind  tunnel  models  could 
be  used  as  tools  to  further  the  imderstanding  of  the 
full-scale  BVI  noise  sources.  Small-scale  wind  tun¬ 
nel  models  provide  a  major  source  of  data  to  validate 
noise  prediction  efforts  and  have  contributed  to  the 
understanding  of  BVI  noise  sources.^”® 

BVI  noise  prediction  techniques  usually  involve  pre¬ 
dicting  the  geometry  of  the  vortices  shed  from  the  heli¬ 
copter  rotor  blade  tips.  This  information  is  then  used 
to  determine  the  blade  loading,  either  by  the  use  of 
lifting  line  theory  or  full-potential  methods.  The  ac¬ 
tual  noise  prediction  methods  are  based  on  the  Ffowcs 
Williams  and  Hawking  equation.®  These  method  re¬ 
quire  that  the  blade  loading  be  supplied  as  input  in 
order  to  predict  the  noise. 

There  is  also  another  method  of  BVI  noise  predic¬ 
tion  being  developed.  This  method  involves  the  use 
of  the  Kirchhoff  equation  to  propagate  the  pressure 
fluctuations  to  the  far-field.  In  this  method,  the  com¬ 
putational  fluid  dynamic  (CFD)  solution  to  either  the 
full-potential  equation  or  the  Euler  equation  is  calcu¬ 
lated  on  the  Kirchhoff  surface.^®  However,  the  rotor 


tip  vortex  geometry  must  be  known  before-hand  as 
input  to  the  CFD  solution. 

Since  all  the  current  prediction  schemes  depend  on 
either  prescribed  or  free- wake  models  of  the  shed  rotor 
tip  vortex  geometry,  and  there  is  little  or  no  data  for 
the  core  size,  vortex  strength  and  location  available 
to  validate  the  wake  models,  it  is  prudent  to  under¬ 
stand  how  the  current  modeling  of  the  tip  vortex  af¬ 
fects  the  BVI  noise  predictions.  One  method  used  to 
evaluate  the  wake  model  is  to  compare  the  locations 
of  the  possible  blade- vortex  interactions  as  predicted 
by  the  wake  model  with  the  locations  of  the  peak  posi¬ 
tive  and  negative  values  in  the  measured  aerodynamic 
surface  pressure.  This  method  determines  how  well 
the  vortex  model  predicts  the  inplane  geometry  of 
the  vortex,  but  does  not  indicate  how  well  the  out 
of  plane  geometry  is  predicted  or  how  well  the  vortex 
core  size  and  strength  are  modeled.  Another  method 
to  evaluate  the  wake  model  is  to  use  the  wake  model 
for  the  prediction  of  the  rotor  aerodynamic  loads  and 
the  acoustics.  The  predicted  aerodynamics  and  acous¬ 
tics  are  then  compared  to  the  measured  data.  In  this 
method,  the  discrepancies  in  the  comparisons  are  due 
to  an  accumulation  of  errors  from  the  vortex  model, 
the  aerodynamic  prediction  and  the  acoustic  predic¬ 
tion.  However,  it  can  be  assumed  from  the  work  doc¬ 
umented  in  Refs.  12-14  that  the  acoustic  prediction 
codes  add  little  error  to  the  end  product.  The  results 
presented  in  these  works  used  measured  aerodynamic 
data  as  input  to  the  acoustic  codes  and  were  in  good 
agreement  with  the  measured  acoustic  data.  There¬ 
fore,  the  comparison  of  the  aerodynamic  predictions 
to  the  measured  data  provides  information  on  how 
well  the  wake  geometry  model  and  the  aerodynamic 
prediction  scheme  are  performing.  The  comparison  of 
the  predicted  acoustics,  using  the  predicted  aerody¬ 
namics,  to  the  measured  acoustics  provides  informa¬ 
tion  about  how  sensitive  the  acoustic  prediction  is  to 
errors  in  the  wake  geometry  and  aerodynamic  models. 

It  is  the  goal  of  the  current  work  to  document  the 
effect  of  the  vortex  modeling  in  both  the  wake  ge¬ 
ometry  prediction  program  and  the  aerodynamic  pre¬ 
diction  program  on  the  BVI  noise  predictions.  Both 
the  wake  geometry  and  the  aerodynamic  prediction 
programs  described  in  this  paper  use  the  Biot-Savart 
Law  to  determine  the  influence  of  the  vortical  veloc¬ 
ities  on  the  blade  loading,  therefore,  the  vortex  core 
size  is  an  input  variable.  A  review  of  previous  work 
documenting  the  effects  of  changing  the  core  size  will 
be  disoissed,  but  the  main  emphasis  will  be  on  the 
presence  of  a  vortex  and  its  location.  Here,  the  vortex 
presence  indicates  whether  or  not  a  vortex  trajectory 
is  included  in  the  aerodynamic  prediction.  Remov¬ 
ing  a  vortex  from  the  prediction  reveals  the  relative 
importance  of  that  vortex  to  the  acoustic  prediction. 

The  wake  geometry  model  used  in  this  paper  is  the 
Scully  model  in  C  AMR  AD/ J  A,  a  comprehensive  an¬ 


alytical  model  of  rotor  aerodynamics  and  dynamics 
by  Johnson  Aeronautics.  FPR,  a  full  potential  ro¬ 
tor  flow  solver,  incorporates  the  wake  geometry  and 
partial  inflow  angles  calculated  by  CAMRAD/JA  in 
the  computation  of  the  aerodynamic  surface  pressures. 
RAPP,  a  rotor  acoustic  prediction  program  based  on 
the  Ffowcs  Williams  and  Hawkings  equation,  calcu¬ 
lates  the  acoustic  pressure  using  the  FPR  computed 
surface  pressure. 

The  rotor  aerodynamic  analysis  in  CAMRAD/JA 
is  based  on  lifting  line  theory  with  separate  models 
for  the  rotor  blade  and  the  wake  that  are  coupled 
by  the  induced  velocity  and  the  bound  circulation. 
The  rotor  blade  model  consists  of  unsteady,  compress¬ 
ible,  viscous  flow  about  an  infinite  aspect-ratio  wing 
and  includes  corrections  for  yawed  flow,  three  dimen¬ 
sional  and  unsteady  effects.  A  lifting-surface  correc¬ 
tion  is  included  for  the  first-order  lifting-line  theory 
blade-vortex  interactions.  The  wake  model  consists 
of  an  incompressible  vortex  wake  from  a  lifting-line 
with  distorted  wake  geometry  (free  wake).  The  in¬ 
board  trailed  and  shed  vorticity  is  modeled  using  vor¬ 
tex  sheets.  The  Biot-Savart  Law  is  used  to  calculate 
the  rotor  wake  induced  velocity.  The  circulation  in  the 
wake  is  determined  by  the  radial  and  azimuthal  vari¬ 
ation  of  the  bound  circulation.^^  The  CAMRAD/JA 
computed  wake  geometry  and  inflow  angles  are  used 
in  FPR  to  solve  for  the  aerodynamic  surface  pressures. 
Since  the  computational  domain  of  FPR  will  include 
some  of  the  wake  geometry,  CAMRAD/JA  calculates 
the  partial  angle-of-attack,  obtained  by  excluding  the 
wake  inside  the  computational  domain  when  evaluat¬ 
ing  the  non-uniform  induced  velocity. 

FPR  solves  the  3-D,  unsteady,  full  potential  equa¬ 
tion  for  transonic  flow  with  the  density  function  deter¬ 
mined  by  the  Bernoulli  equation.  This  rotor  code,  in 
strong  conservation  form,  is  based  on  the  fixed-wing 
code  of  Bridgeman  et  Discrete  vortices  are  in¬ 
troduced  into  the  finite- difference  computations  using 
a  velocity  decomposition  method.  The  total  velocity 
field  is  the  sum  of  the  gradient  of  the  potential  and 
the  known  vortex  velocity.  The  velocities  associ¬ 
ated  with  the  discrete  vortices  are  determined  using 
the  Biot-Savart  Law  as  described  by  Scully  in  Ref.  20. 
The  vortex  is  defined  as  a  series  of  straight  line  seg¬ 
ments  with  the  strength  varying  linearly  along  each 
segment.  An  algebraic  model  is  used  for  the  solid 
body  rotation  of  the  vortex  core.  The  potential  equa¬ 
tion  is  solved  using  half-point  differencing  formulas. 
First-order  backward  differencing  is  used  in  time  and 
second-order  differencing  is  used  in  space.  The  grids 
for  the  solution  are  a  spanwise  series  of  body  fitted 
O-grids  with  80  points  in  the  chordwise  direction,  25 
points  in  the  spanwise  direction,  and  25  points  in  the 
normal  direction.  The  time  step  for  the  solution  cor¬ 
responds  to  0,125°  rotor  azimuth,  or  about  64  time 
steps  per  chord  at  the  rotor  tip.  The  reason  for  this 
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relatively  small  time  step  is  to  capture  the  BVI  events 
that  happen  in  a  very  short  period  of  time. 

A  simple,  yet  accurate  rotor  acoustic  prediction  pro¬ 
gram  (RAPP)  utilizes  the  Ffowcs  Williams  and  Hawk¬ 
ings  (FW-H)  equation  in  a  form  well  suited  to  in¬ 
corporate  blade  surface  pressure  from  computational 
codes  such  as  FPR.  The  Lighthill  stress  tensor,  a  vol- 
xime  integral  that  models  the  noise  produced  by  a 
non-linear  flow  field,  is  not  included  in  the  compu¬ 
tation  performed  by  RAPP.  It  is  assumed  in  the  cur¬ 
rent  work  that  the  blade-vortex  interactions  do  not 
result  in  strong  shocks  and  that  the  contribution  to 
noise  from  this  non-linear  flow  field  term  is  negligible. 
RAPP  can  model  the  force  terms  in  the  FW-H  equa¬ 
tion  as  either  compact  or  noncompact  sources.  When 
the  chordwise  compact  assumption  is  iised,  the  sources 
are  distributed  spanwise  along  the  quarter  chord  of  the 
acoustic  planform.  The  source  strengths  are  then  set 
equal  to  the  aerodynamic  sectional  coefficients  calcu¬ 
lated  by  FPR.  When  the  noncompact  assumption  is 
used,  the  sources  are  distributed  along  the  chord  at 
each  spanwise  station.  The  source  strengths  are  equal 
to  the  local  surface  pressure  calculated  by  FPR.  When 
experimentally  measured  aerodynamic  surface  pres¬ 
sure  is  used  as  input  for  both  the  compact  and  non¬ 
compact  formulations,  there  is  little  difference  in  the 
predictions  of  the  acoustic  pressure.  However,  when 
the  aerodynamic  pressure  from  FPR  is  used  as  input, 
there  is  a  noticeable  difference  between  the  acoustic 
pressure  predicted  by  the  two  formulations.  This  dif¬ 
ference  in  the  acoustic  predictions  will  be  shown  for 
both  the  2-blade,  teetering  hub,  AH-l/OLS  model 
rotor  and  for  the  4-bladed,  fully  articulated,  UH60 
model  rotor. 

Core  Size  Effects 

Since  both  C  AMR  AD/ J  A  and  FPR  use  the  Biot- 
Savart  Law  to  determine  the  influence  of  the  vorti¬ 
cal  velocities  on  the  blade  loading,  both  codes  use 
an  algebraic  vortex  model  for  the  solid  body  rota¬ 
tion  of  the  vortex  core.  Previous  work  documented 
the  changes  in  the  wake  geometry  and  the  partial  in¬ 
flow  angles  predicted  by  CAMRAD/JA  when  the  core 
radius,  normalized  by  the  blade  chord,  was  increased 
from  ajc  =  0.046  to  a/c  =  0.20.  The  work  was 
performed  for  the  AH-l/OLS  model  rotor  system  op¬ 
erating  at  a  flight  condition  with  an  advance  ratio 
equal  to  0.164,  shaft  tilt  angle  equal  to  1.0  degree, 
hover  tip  Mach  number  of  0.665,  and  a  thrust  coeffi¬ 
cient  of  0.00535.  The  in-plane  wake-geometry  did  not 
vary  with  core  size  and  the  out  of  plane  wake  geome¬ 
try  and  vortex  strength  varied  only  slightly  with  core 
size.  The  partial  inflow  angles  became  a  smoother 
function  of  azimuth  for  the  increased  core  size.  In  an 
attempt  to  improve  the  modeling  of  the  wake  geome¬ 
try  by  CAMRAD/JA  at  the  inboard  radial  stations. 


Azimuth  (deg) 

Fig.  1  CAMRAD/JA  partial  angles  as  a  function  of 
azimuth 

the  core  size  was  increased  drastically  to  a/c  =  2.0  in 
the  current  work  as  recommended  by  Johnson  Aero¬ 
nautics.  Even  with  the  large  increase  of  core  size,  the 
in-plane  wake  geometry  did  not  change,  and  the  out- 
of-plane  wake  geometry  and  vortex  strength  changed 
only  slightly.  However,  there  was  a  significant  differ¬ 
ence  in  the  partial  inflow  angles  predicted  by  CAM¬ 
RAD/JA  at  the  inboard  radial  stations.  The  differ¬ 
ence  was  not  as  large  at  the  outboard  radial  stations. 
The  partial  angles  are  plotted  in  Fig.  1  as  a  func¬ 
tion  of  azimuth  for  a/c  =  0.046,0.20,  and  2.0.  Al¬ 
though  it  was  not  realized  in  the  work  of  Ref.  21,  the 
partial  angles  should  not  change  with  core  size.  The 
calculation  of  the  partial  angle  should  remove  any  vor¬ 
tex  influence  from  the  inflow  angle.  This  obviously  is 
not  being  accomplished  and  should  be  noted  by  any 
user  of  CAMRAD/JA  calculated  partial  angles.  The 
plots  of  leading  edge  differential  pressure  as  a  func¬ 
tion  of  azimuth  in  Fig.  2  show  the  difference  in  the 
FPR  calculations  when  using  the  wake  geometry  and 
partial  angles  from  CAMRAD/JA  with  core  sizes  of 
a/c  =  0.20  and  a/c  =  2.0.  The  core  size  used  in 
FPR,  a/c  =  0.20,  was  the  same  for  both  cases.  The 
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Fig.  2  Leading  edge  differential  pressure  as  a  function 
of  azimuth  for  the  AH-l/OLS  model  rotor.  Measured 
data  is  compared  to  FPR  calculations  using  partial 
angles  from  CAMRAD/JA  with  core  sizes  of  a/c  = 
0.20  and  2.0. 


Fig.  3  Leading  edge  differential  pressure  as  a  function 
of  azimuth  for  the  AH-l/OLS  model  rotor.  Measured 
data  is  compared  to  FPR  calculations  with  core  sizes 
of  a/c  =  0.20  and  0.50. 


measured  data  for  this  case  is  also  included  in  Fig. 
2.  The  larger  core  size  in  CAMRAD/JA  improves  the 
prediction  inboard  much  more  than  for  the  outboard 
prediction.  All  the  following  aerodynamic  predictions 
by  FPR  discussed  in  the  current  paper  will  use  input 
from  CAMRAD/JA  with  a  core  size  of  a/c  =  2.0.  This 
does  not  define  the  core  size  in  the  FPR  calculations. 
They  are  defined  separately. 

Changing  the  core  size  in  FPR  has  a  noticeable  ef¬ 
fect  on  the  predicted  aerodynamics  and  acoustics.  The 
differential  pressures  near  the  leading  edge,  calculated 
by  FPR  using  a/c  =  0.20  and  0.50,  are  plotted  as  a 
function  of  azimuth  in  Fig.  3.  The  measured  differen¬ 
tial  pressure  is  also  included  in  Fig.  3.  Increasing  the 
core  size  in  FPR  not  only  reduces  the  amplitude,  but 
also  changes  the  character  of  the  leading  edge  pres¬ 
sure.  This  indicates  that  using  a  larger  core  size  in 
FPR  to  compensate  for  errors  in  miss  distance  defined 
by  the  CAMRAD/JA  wake  geometry  will  not  result 
in  a  better  prediction  of  the  aerodynamics. 


The  effect  of  changing  the  core  size  in  FPR  on  the 
acoustic  predictions  is  shown  in  Fig.  4,  where  the 
acoustic  pressure  is  plotted  as  a  function  of  time  for 
one  blade  passage.  The  measured  acoustic  data  and 
the  predicted  acoustic  data  using  the  compact  formu¬ 
lation  of  RAPP  are  shown  for  four  microphone  loca¬ 
tions.  All  the  microphones  are  located  1.72  rotor  di¬ 
ameters  from  the  rotor  hub.  Microphone  2  is  in  the 
plane  of  the  rotor  at  180  degrees  rotor  azimuth.  Zero 
azimuth  is  define  to  be  over  the  tail  of  the  helicopter. 
Microphone  3  is  30  degrees  below  the  rotor  plane  at 
180  degrees  rotor  azimuth.  Microphones  7  and  9  are 
also  30  degrees  below  the  rotor  plane.  Microphone  7 
is  at  150  degrees  rotor  azimuth  and  microphone  9  is 
at  210  degrees  rotor  azimuth.  The  acoustic  pressure  is 
underpredicted  in  the  rotor  plane  and  over  predicted 
out  the  plane  when  the  core  size  is  a/c  =  0.20  in  FPR. 
When  the  core  size  is  a/c  =  0.50  in  FPR  the  acoustic 
pressure  is  underpredicted  for  all  microphones  posi¬ 
tions.  The  acoustic  predictions  using  the  noncompact 
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a)  normalized  time  b)  normalized  time  c)  normalized  time 

Fig.  4  Acoustic  pressure  is  plotted  as  a  function  of  normalized  time  (1  rev  =  1  unit  of  time),  column  a)  the 
measured  acoustic  pressure,  column  b)  the  acoustic  pressure  predicted  using  the  compact  formulation  with  FPR 
surface  pressure  for  a/c  =  0.20,  and  column  c)  the  acoustic  pressure  predicted  using  the  compact  formulation 
with  FPR  surface  pressure  for  a/c  =  0.50. 


0.0  02  0.4  0.0  0.2  0.4  0.0  0.2  0.4 


a)  normalized  time  b)  normalized  time  c)  normalized  time 

Fig.  5  Acoustic  pressure  is  plotted  as  a  function  of  normalized  time  (1  rev  =  1  unit  of  time),  column  a)  the 
measured  acoustic  pressure,  column  b)  the  acoustic  pressure  predicted  using  the  non-compact  formulation  with 
FPR  surface  pressure  for  a/c  =  0.20,  and  column  c)  the  acoustic  pressure  predicted  using  the  non-compact 
formulation  with  FPR  surface  pressure  for  a/c  =  0.50. 
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,“A-,  measured  interaction  location 
-D-,  -A-,  -O- predicted  interaction  location 
—  epicycloid  interaction  location 


Fig.  6  Leading  edge  differential  pressure  as  a  function 
of  azimuth  for  the  AH-1  /OLS  model  rotor.  The  vortex 
locations  of  interest  are  labeled  for  the  measured  data 
and  the  FPR  prediction  with  a/c  =  0.20 

formulation  in  RAPP  are  shown  in  Fig.  5.  The  mea¬ 
sured  data  and  acoustic  predictions  resulting  from  the 
core  size  variation  in  FPR  are  shown  for  the  same 
microphone  positions  of  Fig.  4.  The  acoustic  pre¬ 
dictions  resulting  from  the  noncompact  formulation 
greatly  underpredict  the  in-plane  acoustic  pressure. 
The  out  of  plane  predictions  for  a/c  =  0.20  in  FPR 
do  not  overpredict  the  positive  pulses  as  is  the  case  for 
the  compact  formulation.  The  positive  peak  is  actu¬ 
ally  slightly  imderpredicted  at  microphones  3  and  9. 
Neither  formulation  is  capable  of  predicting  the  second 
positive  peak  seen  at  microphones  2  and  3.  When  the 
core  size  is  increased  to  a/c  =  0.50  in  FPR  the  acous¬ 
tic  pressure  predicted  by  the  noncompact  formulation 
is  greatly  reduced. 

Vortex  Location 

The  vortex  location  defines  the  azimuth  location  of  the 
interaction  and  the  distance  between  the  blade  and  the 
vortex  at  the  time  of  the  interaction.  The  acoustic  pre¬ 
dictions  are  sensitive  to  both  the  inplane  location  of 
the  interaction  and  the  miss  distance.  The  modeling 
of  the  inplane  vortex  geometry  by  CAMRAD/JA  can 
be  evaluated  by  comparing  the  locations  of  the  posi¬ 
tive  and  negative  peaks  in  the  leading  edge  differential 
pressure  from  the  FPR  predictions  to  the  locations 
of  the  positive  and  negative  peaks  in  the  measured 
leading  edge  differential  pressure.  The  negative  peaks 
indicate  the  location  of  the  BVI  on  the  advancing  side 
and  the  positive  peaks  indicate  the  location  of  the  BVI 


Fig.  7  Blade  vortex  interaction  locations  in  the  plane 
of  the  rotor  for  the  AH-l/OLS  model  rotor. 


on  the  retreating  side.  The  leading  edge  pressures  are 
plotted  in  Fig.  6  for  a  radial  station  of  r/R  =  0.846 
with  the  vortex  locations  of  interest  labeled  as  vl,  v2, 
v3  and  v4.  Both  the  predicted  and  measured  pres¬ 
sures  are  shown  in  Fig.  6  for  the  same  2-bladed  rotor 
case  described  above.  The  locus  of  interaction  loca¬ 
tions  for  each  vortex  is  plotted  in  Fig.  7  along  with 
the  calculation  of  all  the  possible  interactions  of  the 
blade  with  the  epicycloid  representing  the  vortex  tra¬ 
jectory.  The  epicycloid  vortex  trajectory  does  not  take 
into  account  any  convection  of  the  vortex  by  induced 
flow  fields.  The  CAMRAD/JA-FPR  predicted  inter¬ 
action  locations  on  the  advancing  side  of  the  rotor  fall 
between  the  epicycloid  and  the  measured  interaction 
locations.  This  indicates  that  CAMRAD/JA  is  mod¬ 
eling  some  convection  from  flow  of  one  vortex  induced 
on  another,  but  that  this  modeling  is  not  adequate. 
The  incorrect  prediction  of  the  inplane  location  of  the 
BVI  could  account  for  why  the  acoustic  prediction  is 
missing  the  second  peak  of  microphone  3  in  both  Figs. 
4  and  5. 

The  prediction  of  the  wake  geometry  becomes  more 
complicated  for  a  four  bladed  rotor  system.  The 
CAMRAD/JA-FPR  prediction  of  the  leading  edge  dif¬ 
ferential  pressure  at  r/R  =  0.92  is  compared  to  the  ex¬ 
perimentally  measured  leading  edge  differential  pres¬ 
sure  at  r/R  =  0.92  in  Fig.  8  for  the  UH-60  model 
rotor  tested  in  the  DNW.®  The  operating  conditions 
for  the  BVI  test  case  shown  in  Fig.  8  are  an  advance 
ratio  of  0.15,  hover  tip  Mach  number  of  0.64,  and 
thrust  coefficient  of  0.0045.  The  positive  and  nega¬ 
tive  peaks  have  been  identified  as  in  Fig.  6.  There 
is  obviously  a  phase  difference  between  the  identified 
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Azimuth  (deg) 

Fig.  8  Leading  edge  differential  pressure  as  a  function  of  azimuth  for  the  UH60  model  rotor.  The  vortex 
locations  of  interest  are  labeled  for  both  the  measured  data  and  the  FPR  prediction  with  a/c  =  0.20 

peaks  in  the  measured  data  and  the  identified  peaks 
in  the  CAMRAD/JA-FPR  prediction. 

The  locus  of  interaction  locations  for  each  vortex  is 
shown  in  Fig.  9  along  with  the  calculated  positions 
of  all  the  possible  interactions  of  the  blade  with  the 
epicycloid  representing  the  vortex  trajectory.  On  the 
advancing  side  of  the  rotor,  the  predicted  BVI  loca¬ 
tions  correspond  to  the  BVI  locations  of  the  epicycloid 
trajectory.  Whereas  the  measured  BVI  locations  oc¬ 
cur  between  or  cross  over  the  BVI  locations  of  the 
epicycloid  trajectory.  This  indicates  that  there  is  a 
large  effect  of  induced  velocity  on  the  location  of  the 
vortex  interaction  that  is  not  being  calculated  cor¬ 
rectly  by  CAMRAD/JA.  The  effect  of  the  induced  ve¬ 
locity  is  much  larger  for  four  the  4-bladed  rotor  than 
for  the  2  bladed  rotor.  This  is  becaiise  there  are  more 
vortices  for  the  4-bladed  rotor  in  a  similar  region  as 
for  the  2-bladed  rotor.  Therefore,  the  vortices  are  in 
closer  proximity  and  exert  a  greater  influence  on  one 
another. 

The  phase  shift  between  the  measured  BVI  loca¬ 
tions  and  the  predicted  BVI  locations  (see  Fig.  8) 
indicates  that  the  out-of-plane  prediction  of  the  vor¬ 
tex  location  by  CAMRAD/JA  is  not  correct.  Fig. 

10  displays  the  predicted  wake  geometry  within  the 
computational  domain  when  the  rotor  is  at  55  degrees 
azimuth.  The  segments  of  the  vortex  trajectories  are 
projected  both  in  the  plane  of  the  rotor  and  in  a  plane 
perpendicular  to  the  rotor.  The  rotor  blade  and  com¬ 
putational  domain  are  also  shown  in  Fig.  10.  The 


-A-,  meassured  interaction  location 

HD-,  -A-,  -0-,  — O—  predicted  interaction  location 
—  epicycloid  interaction  location 

Fig.  9  Blade  vortex  interaction  locations  in  the  plane 
of  the  rotor  for  the  UH60  model  rotor. 
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Fig.  10  C  AMR  AD/ J  A  predicted  wake  geometry  at  55  degrees  azimuth  for  the  UH60  model  rotor. 


Fig.  11  Acoustic  pressure  is  plotted  as  a  function  of  normalized  time  (1  rev  =  1  unit  of  time),  row  a) 
measured  acoustic  data  compared  to  acoustic  prediction  using  compact  and  non-compact  formulations  with 
measured  surface  pressure,  row  b)  measured  acoustic  data  compared  to  acoustic  prediction  using  compact  and 
non-compact  formulations  with  FPR  predicted  surface  pressure. 
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comparison  of  measured  and  calculated  leading  edge 
pressure  (see  Fig.  8)  suggests  that  vortices  v2  and  v3 
in  Fig.  10  are  too  close  to  the  rotor  whereas  vortices 
vll  and  vl2  are  too  far  from  the  rotor. 

The  acoustic  predictions  for  this  UH-60  case  are 
compared  to  the  measured  data  in  Fig.  11  for  mi- 
crophone  9  which  is  3.0  rotor  radii  from  the  hub  and 
25  degrees  below  the  rotor  plane  at  150  degrees  rotor 
azimuth.  Both  the  compact  and  noncompact  formu¬ 
lations  have  been  included  for  the  predictions  using 
measured  surface  pressure  and  for  the  predictions  us¬ 
ing  FPR  calculated  surface  pressure.  Both  the  com¬ 
pact  and  noncompact  formulations,  when  using  mea¬ 
sured  surface  pressure,  compare  well  with  measured 
data.  However,  both  acoustic  formulations  using  the 
FPR  calculated  surface  pressure  overpredict  the  BVI 
noise.  The  poor  comparison  to  the  measured  data  is 
a  direct  result  of  the  incorrect  modeling  of  the  wake 
geometry. 

To  see  how  sensitive  the  acoustic  predictions  are  to 
the  miss  distance  at  the  time  of  interaction,  the  out-of- 
plane  location  predicted  by  C  AMR  AD/ J  A  was  shifted 
down  0.25  chords  and  0.5  chords  then  used  by  FPR 
to  predict  the  surface  pressures.  The  effect  of  miss 
distance  on  the  prediction  of  the  surface  pressures  is 
shown  in  Fig.  12  for  the  2-bladed  AH-l/OLS  rotor 
case  discussed  above  at  radial  stations  of  r/R  =  0.846 
and  0.91  and  a  chordwise  station  of  x/c  =  0.03.  The 
FPR  results  shown  in  Fig.  12  include  the  predictions 
using  the  unshifted  C  AMR  AD/ J  A  wake  geometry  and 
the  predictions  using  the  CAMRAD/JA  wake  geom¬ 
etry  shifted  0.25  chords  and  0.5  chords.  The  vortex 
labeled  vl  in  Fig.  12  is  only  slightly  affected  by  the 
shift  in  the  vortex  geometry  at  either  radial  station. 
Vortices  v2,  v3  and  v4  are  all  significantly  affected  by 
the  geometry  shift.  The  amplitude  of  vortex  v2  is  re¬ 
duced  while  the  amplitude  of  v3  is  increased  because 
v2  was  below  the  plane  of  the  rotor  and  was  moved 
farther  from  the  blade  with  the  downward  shift  and 
v3  was  above  the  plane  of  the  rotor  and  was  moved 
closer  to  the  blade  with  the  downward  shift.  The  effect 
of  changing  the  miss  distance  on  the  predicted  acous¬ 
tic  pressure  is  shown  in  Fig.  13,  where  the  compact 
acoustic  formulation  was  used,  and  in  Fig,  14,  where 
the  non-compact  acoustic  formulation  was  used.  For 
both  formulations,  the  major  acoustic  pulse  is  reduced 
while  the  following  pulses  are  increased.  The  increase 
in  the  secondary  pulses  is  most  evident  for  the  0.25 
chord  shift.  The  predictions  in  Figs.  13  and  14  indi¬ 
cate  that  the  acoustic  predictions  are  quite  sensitive 
to  small  shifts  in  the  wake  geometry. 

Vortex  Presence 

Being  able  to  determine  which  blade-vortex  inter¬ 
action  produces  the  most  noise  is  valuable  knowl¬ 
edge  that  could  lead  to  major  noise  reduction.  One 


Fig.  12  Leading  edge  differential  pressure  as  a  func¬ 
tion  of  azimuth  for  the  AH-l/OLS  model  rotor.  The 
measured  data  is  compared  to  FPR  predictions  for  an 
unshifted  wake  and  for  the  wake  shifted  away  from  the 
rotor  0.25c  and  0.50c. 
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Fig.  13  Acoustic  pressure  is  plotted  as  a  function  of  normalized  time  (1  rev  =  1  unit  of  time),  column  a)  the 
measured  acoustic  pressure,  column  b)  the  acoustic  pressure  predicted  using  the  compact  formulation  with  FPR 
surface  pressure  with  the  wake  shifted  0,25c,  and  column  c)  the  acoustic  pressure  predicted  using  the  compact 
formulation  with  FPR  surface  pressure  with  the  wake  shifted  0.50c. 
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Fig.  14  Acoustic  pressure  is  plotted  as  a  function  of  normalized  time  (1  rev  =  1  unit  of  time),  column  a) 
the  measured  acoustic  pressure,  colunm  b)  the  acoustic  pressure  predicted  using  the  non-compact  formulation 
with  FPR  surface  pressure  with  the  wake  shifted  0.25c,  and  column  c)  the  acoustic  pressure  predicted  using  the 
non-compact  formulation  with  FPR  surface  pressure  with  the  wake  shifted  0,50c. 


X  (aspect  ratio) 


y  (aspect  ratio)  radius  (aspect  ratio) 

Fig.  15  CAMRAD/JA  predicted  wake  geometry  at  70  degrees  azimuth  for  the  AH-l/OLS  model  rotor. 

technique  that  is  available  only  to  the  prediction  Removing  the  second  vortex  trajectory  removes  the 
method  is  to  remove  vortices  from  the  computation.  largest  interaction  from  the  FPR  predictions  shown  in 
The  extent  that  the  acoustic  predictions  are  affected  Fig.  16  and  also  removes  the  largest  interaction  from 

by  the  vortex  removal  indicates  the  relative  impor-  the  acoustic  predictions  in  Figs.  17  and  18.  This  ex- 

tance  of  that  vortex  to  the  noise  signature.  This  proce-  cersize  illustrates  the  ability  to  identify  the  BVI  noise 

dure  was  preformed  for  the  2-bladed  AH-l/OLS  model  sources, 
rotor  discussed  previously.  There  are  three  vortex 
trajectories  included  in  the  computation  performed 

by  FPR.  The  segments  of  these  trajectories  that  are  Concluding  Remarks 

within  the  computational  domain  when  the  rotor  is 

at  70  degrees  rotor  azimuth  are  shown  in  Fig.  15,  Most  of  the  current  BVI  noise  prediction  schemes  rely 

which  is  a  snapshot  of  the  vortex  segments  at  a  single  on  prescribed- wake  or  free- wake  models  to  define  the 

instant.  The  purpose  of  Fig.  15  is  to  clarify  which  vortex  location  and  strength  of  the  tip  vortices.  The 

vortex  is  being  removed  from  the  computation.  Fig.  current  work  shows  that  there  can  be  no  compensating 

16  shows  the  effect  of  removing  the  vortices  from  the  for  the  discrepancies  in  the  wake  model  when  attempt- 

computation  on  the  leading  edge  differential  pressure.  iug  to  accurately  predict  the  BVI  noise.  More  specif- 

The  vortices  have  been  labeled  to  match  the  label-  ically,  errors  in  modeling  the  vortex  location  can  not 

ing  of  Fig.  6.  The  interactions  can  be  removed  from  be  compensated  for  by  adjusting  the  core  size  in  the 

the  calculated  pressure  by  removing  a  vortex  trajec-  viscous  core  model.  This  paper  discusses  the  effects 

tory  from  the  computation.  In  two  cases,  removing  of  modeling  the  core  size  in  both  CAMRAD/JA  and 

one  vortex  trajectory  actually  removes  an  interaction  FPR  on  the  acoustic  predictions.  The  most  important 

from  both  the  advancing  side  and  the  retreating  side,  result  of  varying  the  core  size  in  CAMRAD/JA  is  that 

indicating  that  the  vortex  trajectory  has  segments  in  CAMRAD/JA  does  not  accurately  remove  the  vortex 

the  computational  domain  on  both  sides  of  the  rotor  influence  from  the  inflow  angles  when  calculating  the 

partial  angles.  Otherwise,  changing  the  core  size  in 
CAMRAD/JA  has  little  affect  on  the  wake  geometry 
The  effect  of  removing  the  vortex  trajectories  from  or  the  predicted  acoustics.  The  aerodynamic  surface 

the  computation  on  the  acoustic  predictions  is  shown  pressure  predicted  by  FPR  is  quite  sensitive  to  the 

in  Figs.  17  and  18.  In  Fig.  17,  the  compact  for-  core  size  used  for  the  predictions.  This  in  turn  greatly 

mulation  is  used  to  predict  the  acoustics  and  in  Fig.  affects  the  acoustic  predictions  that  use  the  FPR  pre- 

18,  the  non-compact  formulation  is  used.  These  plots  dieted  surface  pressure.  For  both  the  2-bladed  and 

can  be  compared  to  the  first  two  colunms  of  Figs.  4  4-bladed  rotor  systems  discussed  in  this  paper,  CAM- 

and  5  which  show  the  measured  acoustic  data  and  the  RAD/JA  did  not  adequately  model  the  convection  of 

acoustic  predictions  using  the  same  core  size  and  vor-  the  vortex  trajectories  caused  by  induced  flow  that  is 

tex  positions  used  for  the  vortex  removal  calculations.  produced  by  nearby  vortex  trajectories.  The  error  in 
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Fig.  16  Leading  edge  differential  presssure  as  a  function  of  time  for  the  AH-l/OLS  model  rotor.  The  FPR 
predictions  are  shown  with  all  the  vortices  present  in  the  computation  and  with  vl  removed,  v2  and  v5  removed, 
and  v3  and  v4  removed. 


modeling  the  vortex  convection  produces  errors  in  the 
vortex  location  with  respect  to  the  blade.  The  results 
of  this  paper  show  that  both  the  aerodynamic  and 
acoustic  predictions  are  quite  sensitive  to  small  errors 
in  the  vortex  locations. 

The  BVI  noise  prediction  schemes  can  be  quite  use¬ 
ful  when  they  work  relatively  well.  The  current  work 
shows  that  the  inclusion  or  exclusion  of  a  vortex  in 
the  FPR-RAPP  calculation  allows  for  the  determina¬ 
tion  of  the  relative  importance  of  that  vortex  as  a 
BVI  noise  source.  Being  able  to  identify  which  vortex 
is  responsible  for  the  largest  noise  source  will  lead  to 
improved  BVI  noise  reduction  methods. 
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a)  nonnalized  time  b)  normalized  time  c)  normalized  time 


Fig.  17  Acoustic  pressure  predicted  using  the  compact  formulation  is  plotted  as  a  function  of  normalized  time 
(1  rev  =  1  unit  of  time),  column  a)  vl  removed  from  the  FPR  computation, column  b)  v2  and  v5  removed  from 
the  FPR  computation,  and  column  c)  v3  and  v4  removed  from  the  FPR  computation. 


a)  normalized  time  b)  normalized  time  c)  normalized  time 

Fig.  18  Acoustic  pressure  predicted  using  the  noncompact  formulation  is  plotted  as  a  function  of  normalized 
time  (1  rev  =  1  unit  of  time),  column  a)  vl  removed  from  the  FPR  computation,  column  b)  v2  and  v5  removed 
from  the  FPR  computation,  and  column  c)  v3  and  v4  removed  from  FPR  computation. 
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